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Abstract 

Context. One of the key ingredients in establishing the relation between input signal and output flux from a spectrometer is 
accurate determination of the spectrophotometry calibration. In the case of spectrometers onboard satellites, the accuracy of 
this part of the calibration pedigree is ultimately linked to the accuracy of the set of reference spectral energy distributions 
(SEDs) that the spectrophotometric calibration is built on. 

Aims. In this paper, we deal with the spectrophotometric calibration of infrared (IR) spectrometers onboard satellites in the 2 
to 200 fim wavelength range. We aim at comparing the different reference SEDs used for the IR spectrophotometric calibration. 
The emphasis is on the reference SEDs of stellar standards with spectral type later than AO, with special focus on the theoretical 
model atmosphere spectra. 

Methods. Using the MARCS model atmosphere code, spectral reference SEDs were constructed for a set of IR stellar standards 
, (A dwarfs, solar analogs, G9-M0 giants). A detailed error analysis was performed to estimate proper uncertainties on the 
predicted flux values. 

Results. It is shown that the uncertainty on the predicted fluxes can be as high as 10%, but in case high-resolution observational 
optical or near-IR data are available, and IR excess can be excluded, the uncertainty on medium-resolution SEDs can be 
reduced to 1-2% in the near-IR, to ~3% in the mid-IR, and to ~5% in the far-IR. Moreover, it is argued that theoretical 
stellar atmosphere spectra are at the moment the best representations for the IR fluxes of cool stellar standards. 
Conclusions. When aiming at a determination of the spectrophotometric calibration of IR spectrometers better than 3%, effort 
should be put into constructing an appropriate set of stellar reference SEDs based on theoretical atmosphere spectra for some 
, 15 standard stars with spectral types between AO V and M0 III. 

Key words. Instrumentation: spectrographs, Techniques: spectroscopic, Stars: atmospheres, Stars: late-type, Infrared: stars 



1. Introduction 

For (IR) spectrometers, three main steps have to be cov- 
ered to establish the relation between input signal, flux, 
and wavelength: spatial, spectral, and photometric cali- 
bration. In this paper, we focus on the reference spectral 
energy distributions (SEDs) used in the spectrophotomet- 
ric ( — spectral+photomctric) calibration process, with 
main emphasis on the spectral calibration part. 

For (IR) spectrometers onboard satellites, the determi- 
nation of the relative spectral response function (RSRF), 
which characterises the wavelength-dependent response of 
a spectrometer, is often a two-step process: (1) first, the 
RSRF is determined during laboratory tests from mea- 
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suring a cryogenic blackbody calibration source at differ- 
ing temperatures; (2) after launch, the RSRF is refined 
by comparing observations with reference SEDs of vari- 
ous celestial calibration sources, preferably at a spectral 
resolution comparable to that of the instrument. 

Whenever an instrument covers a new spectral window 
or has a higher sensitivity than preceding instruments, 
new reference SEDs have to be constructed. This was and 
will be the case for instruments onboard the ESA- Infrared 
Space Observatory (ISO), the NASA-Spitzer satellite, 
the ESA-Herschel mission, the NASA- James Webb Space 
Telescope (JWST), and many others. In this paper, em- 
phasis will be put on the reference SEDs of calibration 
sources for instruments covering the 2 to 200 /im wave- 
length range. When selecting calibration sources, issues 
like brightness in the wavelength regime covered by the 
instrument, confusion by neighbouring sources, sky visi- 
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bility, and pointing accuracy are to be considered. In the 
2-200 fim wavelength range, planets, asteroids, and cool 
standard stars are the ideal calibration sources to cover 
the whole dynamic range of present-day instruments (see 
Sect. [2]). This paper focuses on the reference SEDs of stan- 
dard stellar candles, with the aim of assessing the reliabil- 
ity and achieving accuracy when using IR stellar standards 
for spectral calibration purposes. Quite often, calibration 
scientists or observers (have to) rely on 'quoted' errors, 
without having insight into this part of the calibration 
pedigree. Since for a proper interpretation of the output 
signal one needs to have a grip on the uncertainties that 
propagate to the fluxes ultimately calculated by an instru- 
ment, the discussion in this paper is required. 

Section [3] gives an overview of the different IR flux 
calibration sources. The status of stellar reference SEDs 
in the IR is discussed in Sect. [3J and we focus on the 
theoretical atmosphere spectra in Sect. [4] Different stellar 
reference SEDs are compared in Sect. [5l and we end with 
the conclusions in Sect. [71 



2. Infrared flux calibration sources 

In this section, we discuss different IR calibration sources 
used in the spectrophotometric calibration process. 



2.1. Planets 

For calibration purposes at wavelengths longer than 
~30 /mi, planets are extremely useful since they are among 
the few astronomical objects bright enough in the far-IR 
to allow sufficiently accurate flux density predictions. At 
100 /im, the flux values for Neptune and Uranus range 
between 100 and 1000 Jy. Uranus is often used as the pri- 
mary calibrator, as it i s known to be a re liable calibrator 
((Griffin fe OrtorJll993tlSidher et al.ll2003h . while Neptune 
and Mars are used as secondary calibrators. 

" 19981: 



Good models for Uranus exist 



Orton fc BurgdorJ l2003t ISidher"et all ^00^^ Radiative 



Moreno 



transfer calculations are done in spherical geometry and 
take the limb into account. The models do not, how- 
ever, take chemical reactions into account, nor heating 
and cooling. In the mm and submm-regimes, the mod- 
els are confirmed within the accuracy of 10 - 20%. The 
big advantage of Mars over Uranus is that it is bright 
Different models for Mars' surface exist 
200fl 



1987; Lellouch et al 



Hart ogh et al 



. .Rudy et al, 
2005J). Using 



Mars as a primary calibrator, however, would pose dif- 
ferent problems: (1) Mars is a planet with a surface and 
an atmosphere, both contributing to the continuum emis- 
sion, (2) sandstorms influence the line shapes, (3) the ice 
caps can influence the continuum if the ice is melting. 
Neptune could serve almost equally well as an independent 
(secondary) calibrator. However, many large features are 
visible in the Voyager IRIS spectra, which may indicate 



2.2. Asteroids 

Typical surface temperatures of main-belt asteroids are 
such that they emit the bulk of their thermal radiation 
in the far-IR. At wavelengths longer than ~20/zm, the 
largest asteroids are brighter than the brightest IR, stellar 
sources: they cover the flux range between 100 and 1000 Jy 
in the wavelength range between 30 and 45 /im, and hence 
fill a gap where stellar calibrators are not available (see 

Fig. 1 in iMiiller fc Lagerrosl[T998I). 

T he work of Muller feller fc Lagerrosl Il998l 12002 . 

20031 : iMuller et all l2005h has largely improved the ac- 
curacy of the theoretical model SEDs of asteroids. The 
most accurate SEDs at present are based on the thermo- 
physical mod el (TPM) for describing t he asteroids' ther- 
mal emission ( Muller fc Lagerrosl 1998f l. Due to the astcr- 
oids' orbit, brightness te mperatures vary by 17 ± 13%. 
Muller fc Lagerrosl (|2002h quote an accuracy of 5 % be- 
tween 5 and 200 /jm for Ceres, Pallas, and Vesta, i.e. for 
the three largest asteroids. 

The main limitations for asteroids as IR calibration 
sources come from the changing background conditions 
and the flux changes on timescales of hours due to rota- 
tion. Moreover, the list of reliable asteroids is still short. 



2.3. Stars 

In the ncar-IR, stars are ideal calibrators, especially since 
they are almost point-like and span a range in flux level 
more than 4 orders of magnitude between 2 and 50 /im, 
and one can create a database covering the whole sky. 
Stellar standards are, however, quite faint in the far-IR 
and can henceforth not be used as primary calibrators at 
the far-IR wavelength ranges. Good IR calibration sources 
need to comply several criteria: be single and non- variable, 
not have an IR excess due to a chromosphere, debris disk, 
or a circumstellar envelope, and be located in an un- 
crowd ed region that ca n be observed all or most of the 



time (jDecin et al.l 120071 ). Preferably, the IR standard is 



syste matic variability with rotational phase ((Bishop et al 
19981) . 



cooler than ~ 10 000 K to provide a good signal-to-noise 
ratio over most of the wavelength range. 

In the case of spectrophotometric calibration, one ex- 
tra requirement should be added to this list regarding 
the composition of the sample of used calibrators with 
regard to the spectral types. For photometric calibration 
purpos es, the use of on ly one spectral type can be jus- 
tified ( Diaz-Miller 20071 ). but in the case of a spectro- 
metric calibration, where every single wavelength needs 
to be calibrated, none of the available reference SEDs of 
any spectral type has high enough accuracy in the full 2 
to 200 /xm wavelength range (see Sect[3|). Three classes of 
spectral stellar standards have been commonly used in the 
IR spectrometric calibration pedigree: (i) early A dwarfs, 
(ii) solar analogs, and (Hi) late- type giants, usually of 
spectral type G9 - M0 III. Each of these groups provides 
a different challenge (see below), but by combining them, 
they (i) reduce the chances for systematic errors, possibly 
introduced by the use of only one spectral type, and (ii) 
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they will increase the total accuracy since the best SED 
can be used for each part of the spectrum. 

3. Status of stellar reference SEDs in the IR 

3.1. Black body 

Planck curves represent the simplest way to model stellar 
far-IR fluxes. A comparison with the predictions of a solar 
continuum model shows clear deviati ons, arising to more 
than 20% in the far-IR (see Fig. 14 in Ivan der Bliek et al. 
19961) . The reason is the contribution of the H free-free 
opacity, shifting the flux-forming region more to outer 
cooler layers. Neither line or continuous absorption nor 
back- warming effects arc/can be included. 

3.2. Engelke function 

The Engelke function (jEnge is a more sophisti- 



cated two-parameter analytical approximation to the 2 - 
60 yum infrared continuum spectrum for giants and dwarfs 
with effective temperature, T e g, between 3500 and 6000 K. 
This spectral function is based on the scaling of a semi- 
empirical (plane-parallel) solar atmospheric profile to dif- 
fering effective temperatures. The result describes the con- 
tinuum spectrum expected for stars in terms of their effec- 
tive temperature and angular size. The estimated probable 
error in absolute flux values is quoted to be ~3 % below 
10 /mi, growing to ~5 % around 25 /tm and 6 % at 60 /zm. 

The main limitation of the Engelke function is neglect 
of the influence of the surface gravity. In Fig. Q] the dif- 
ference between the Engelke function and the continuum 
flux of a plane-parallel theoretical stellar atmosphere (see 
Sect. 2]) with an effective temperature, T e ff, of 3500 K for 
different values of the gravity is shown. For values of the 
logarithm of the gravity, log g, lower than 3.00, the abso- 
lute deviations may be as high as 25 %. At a temperature 
of 6000 K, the absol u te dev iations approximate the errors 
quoted by lEneelkd dl992h . Important, however, for the 
discussion in this paper of the RSRF determination are 
the systematic differences between the Engelke continuum 
and both (1) the continua calculated by ab-initio model 
atmosphere calculations (see Fi g, ffl) and (2) real co ntinua 



Engelke et aT1l2006h . This 



deduced from observations (see . 
systematic discrepancy occurs both for solar analogs and 
the cooler K-M giants: the same type of curves as dis- 
played for the lower gravity values in Fig. Q] (M-giants) 
also occur in an analogous plot for G dwarfs, but with a 
smaller amplitude of ~ 3 %. Note also that for almost all 
values of log g, the ratio as displayed in Fig. [T] is system- 
atically higher than one, implying that when either the 
effective temperature or the angular diameter (8d) is kept 
fixed, the other parameter will be overestimated. 

Another issue not taken into account are sphericity 
effects, important when dealing with the extended atmo- 
spheres of giants. Figure [2] compares a theoretical atmo- 
sphere spectrum (see Sect. H|) calculated using spherical 
geometry with the spectrum obtained in plane-parallel ge- 




Wavelength [micron] 

Figure 1. Ratio between the continuum flux predicted 
from a theoretical plane-parallel stellar atmosphere with 
T e ff = 3500 K and values for the logarithm of the gravity 
ranging between 5.00 and 0.00 and the flux computed from 
the Engelke function. 



ometry for the case of an M0 III giant. Compared to a 
plane-parallel atmosphere, the radiation field of a spheri- 
cal model is diluted in the upper photospheric layers, caus- 
ing the temperature (and hence the source function) there 
to be lower. The lower surface flux is, however, compen- 
sated by the fact that the IR flux arises from higher lay- 
ers. The net result is an infrared excess of the spherical 
model relative to the plane-parallel model. Important for 
this paper is the rising continuum in the bottom panel 
of Fig. [2j implying that the use of the Engelke function 
or a plane-parallel geometry to represent the SED of gi- 
ants may introduce uncertainties on the broad-band RSRF 
characterisation of a few percent. 
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Figure 2. Upper panel: Model atmosphere spectrum in 
spherical geometry (black) compared to the theoretical 
spectrum calculated in plane-parallel geometry (grey) 
representing the M0 III giant (5 And (A/AA = 100). 
Bottom panel: Ratio between spherical and p lane-parallel 
spectru m. Stellar parameters are taken from iDecin et al.l 
(l2003ah . 
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3.3. The IR calibration framework by Cohen et al. 

In a series of papers on spectral irradiance calibration in 
the infrared, Cohen and colleagues have presented a self- 
consistent network of absolutely calibrated reference spec- 
tra between 1.2 and 35 fim for over 600 stars spanning 
spectral types from AO to M 4 and distributed over the en- 
tire sky (ICohen et alll992al lbl. fl995Lll996aL[l99flLl200^L of- 
ten referred to as the CWW network) . The CWW network 
is the result of a tremendous and consistent work done 
during the past 15 years. Four main steps can be distin- 
guished in the network, (i) The approach is based on a pair 
of absolutely calibrated mod els of the two A - type dwarfs 
Vega a nd Sirius calculated bv lKurucd(|l993l) ( Co hen et all 
1992al) . (ii) The next layer in the network is a set of 



composites of 13 secondary calibration stars. Composites 
are spectra measured by various ground-based and air- 
borne telescopes (KAO, IRAS LRS, etc.) that are averaged 
and/or spliced together to form a continuous spectrum be- 
tween 1.2 and 35 fim. This set compromises one G2 dwarf 

- q 1 Cen (G2 V); seven K giants — a Tau (K5 III), a 
Boo (K2 IIIp), 13 Gem(K0 III), a Hya (K3 II-III), a TrA 
(K2 III), e Car (K3 III), 7 Dra (K5 III); and 5 M giants 

- f3 Peg (M2.5 II-III), (3 And (MO II I), a Cet (Ml 5 III) 



7 Cru (M3. 4 III), n UMa (MO III) (jCohen et alj|l992b . 
1995 lll996albh . The KAO was no longer operational at the 
time the composites of 7 Dra, a Cet, 7 Cru and \i UMa 
were constructed. Spectral regions that are opaque from 
the ground were therefore replaced by spectra of other 
stars with the same spectral type. Note that the 'compos- 
ite' of a 1 Cen is a Kurucz model-atmosphere spectrum, 
which was then used as a reference to create the composite 
spectra of a TrA and e Car. For wavelengths longer than 
~22 /im, the Engelkc function was used for twelve compos- 
ites to represent the stellar spectrum, the exception being 
a Tau. (iii) These 13 composite spectra are the basis for a 
set of absolutely calibrated spectral templates for 602 stars 



with spect ral types between G9.5-M0.5 III (| Cohen et al 
1999ll2003h . For each star a smoothed composite spectrum 



is chosen as its spectral template according to the spectral 
class. The template is corrected for reddening and nor- 
malised using available optical photometry. For purpose 
of the NASA-Spitzer Infrared Array Camera (IRAC) cali- 
bration an extra set of absolutely calibrated 0.275 — 35 /im 
spectra of 33 optical standard stars with spectral types be- 



tween AO 
to V ~ 11 



A5 V and K0 - MO III, 



and fl ux leve ls down 
120031) . 



12, was prepared (|Cohen et al.ll2003| ). (iv) In 
support of the Infrared Space Observatory (ISO) PHT and 
LWS instruments, CWW have provided extr apolated con 



tinuum spectra of few composites till 300 /im (j Cohen et al 
1996b| l. The continuum spectra have been interpolated 



in a grid of four plane-parallel model-atmosphere spec- 
tra (see also Sect. I3.2[) with solar abundances. To obtain 
the temperature-versus-continuum optical depth relation 
(r cont ), a scaling based on the ratio T e g (desired) /T e g (grid 
model) was performed. Differences in gravity were not ac- 
counted for. An uncertainty of 5.67 % was assessed to these 
extrapolations. 



It should be emphasised that the CWW network 
has some major advantages, the main ones being (1) 
the common calibrati on pedigree, and (2) the good sky- 
coverage (see Fig. 2 in Cohen et al.ll 1999^ ) . For that reason, 
the CWW network was used to support the calibration 
of many spacebornc, airborne, and ground-based instru- 
ments, as e.g. instruments onboard ISO, Spitzcr, AKARI, 
MSX, etc. 

However, aiming for a (spectrophotometric) calibra- 
tion better than a few percent, some steps used in the con- 
struction of the CWW network deserve some comments. 



A first issue concerns the danger of propagating of 
systematic errors from the top of the calibration 
pedigree. In this context, we want to note that the 
adopted spectra for the fundamental reference stars 
in the calibration pedigree are absolutely-calibrated 
Kurucz models of the A-type stars Sirius and Vega. 
Since no absolute flux measu rements were available 
for Sirius, Cohen et al.l (1992a) chose to calibrate the 
Sirius model absolutely with respect to Vega by us- 
ing the observed magnitude differences in K, L, L', 
M, and at four mid-IR wavelengths. However, the use 
of Vega (or its model spectrum representing an 'ideal' 
Vega without flux ex cess) is stro ngly debated (for an 
overview, we refer to lGravll2007l ). We summarise that 
(i) it is obvious that the Vega debris disk system con- 
trib utes to the flux excess fro m this star beyond 10 fim 
(e.g.lNeugebauer et al.lll984h , (ii) that recent results by 
Absil et all (|2006l ) show an excess of 1.29 ± 0.19% in 
the K-band close-in to Vega (within 1") presumably 
due to a hot inner circumstellar disk; (iii) and Vega 
is a rapidly rotating polc-on star with a significant 
temperature gr adient (~1500 - 2250 K) from its pole 



to its equator (jGulliver et al.l Il994l ; lAufdenberg et al 
2006). so its SED differs from conventional models. 



not assum i ng a r otating star (|Gravll2OO70 . As noted 
by iBohlinl (|2007l) a single-temperature T eff = 9400 K 



Kurucz model fits the Hubble Space Telescope (HST) 
STIS SED to 1-2% from 3200 to 10000 A. However, a 
multi-temperature model may be required to achieve 
a ~1 % precision for shorter and longer wavelengths 
(the topic of this study). Finally (iv) previous work 
has suggested an inconsistency between infrared and 



dard AO V star model (e.g. Rieke et al. 


19851 Megessier 


1995; Verhoelstll2005: Rieke et al.ll2007 


) , probably hav- 



ing its roots in Vega being a rapidly rotating star seen 
pole-on. When the absolute calibration of the Sirius 
spectrum by magnitude-difference with Vega is done 
in the optical and not in the IR, t he fluxes of the 
Sirius spectrum i ncreases by 



3.7% (IVerhoelstl 12005 
Rieke et al. 20071) c ompared to the abs olute calibra- 
tion as proposed bv lCohen et al. 1 1992a ) , shifting the 



angula r diameter of Sirius from 6.04 mas (jCohen et al 



1992al) to 6.15 mas (see Fig. [3]). Based on the Midc ourse 
Space Experiment (MSX) data. lPrice et alj (2004) pro- 
pose an increase in the SED of Sirius by 1 %. This 
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increase in the flux density of the primary standard 
Sirius has a direct influence on the absolute flux level 
of all composites and templates, whose absolute photo- 
metric data are based on relative flux differences com- 
pared to the primary standards. 



Wavelength (^m) 

Figure 3. Observed (red - diamonds) versus synthetic 
(green - stars) magnitude difference between Vega and 
Sirius (Magvcga — Magsh-ius) in the Geneva ph otomet- 
ric system and the JHKL bands ( Verhoelstl 2005 ). Black 
dotted line represents t he difference betwe en both of 
Kurucz's models used in ICohen et all (|f992ah . with stel- 
lar parameters T cff = 9550K, log g = 3.95, [Fe/H] = -0.3 
for Vega and T cff = 9850K, log g = 4.25, [Fc/H] =+0.5 for 
Sirius. However, the Kurucz's model for Sirius is scaled 
in such a way as to match the observed difference in the 
Geneva photometric system, yielding an angular diame- 
ter for Sirius of 6.1 5 mas (instead of 6.04 mas as used in 
Cohen et al1ll992al) . Clearly, the (N)IR magnitude differ- 
ences are not compatible asthey differ by up 10%. 

— One should be aware that, for all composite spectra, 
with the exception of a Tau, the flux values between 
~ 22 and 35 fim are represented by the Engelke func- 
tion. As discussed in the previous section, the neglect 
of both gravity and sphericity effects results in a sys- 
tematic broad-band inaccuracy of ~ 3 % for early A 
dwarfs, risin g to ~ 8 % for late K-giants. ICohen et al 



( 1996bl ) and ICohenl (|2007t ) justify this assumption by 
comparing the UKIRT CGS3 16 - 24 [im spectra of 
a Boo and a Cet with their relevant Engel ke approx- 



imation (see Fig. 5 in Cohen et al. 1996bl) and f rom 



the KAO spectrum of a Tau (|Cohen et al.|[l992bh . A 
comparison to these observations in the A 4 F\-space 
(Fig. [5]) does, however, not give any more support to 
the "flatness" of an Engelke spectrum than to the ris- 
ing continuum of a theoretical model atmosphere (see 
also Sect. [5]). One can only argue that, from a phys- 
ical point of view, the increase in the H opacity to- 
wards longer wavelengths (cr(Hg) oc A 2 ) yields a shift 
of the flux-forming region towards higher, cooler, at- 
mospheric layers, and hence a decrease in the stellar 
surface flux. 
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Figure 4. UKIRT CGS3 15 - 22 /jm spectrum for a Boo 
(K2 IIIp; unsealed, dash-dot line) and KAO spectrum for 
a Tau (K5 III; scaled by a factor 1.11, dashed line) are 
compared with the Enge lke approxim a tion (fu ll line) with 
parameters as used by ICohen et al.1 (|l996bh for a Boo 
(T ff = 4362K, 6d = 21.12mas), and an appropriate the- 
oretical MARCS model at a resolution A/ A A = 1500 (grey 
line) and at a resolution of A/ A A = 70 (dotted line; cor- 
responding to the UKIRT C GS3 resolution), wit h stellar 



parameters as determined in [Dec in et al. | (l2003al) . 



The assumption in the creation of a 'template' spec- 
trum is that a de-reddened spectrum of a star is repre- 
sentative of the infrared spectral signature of any star 
of the same two-dimensional spectral classification. 
This is only true to some extent. Differences in stellar 
parameters and photosphcric abundances can result in 
substantial differences in the strength of molecular ab- 
sorption bands as can be seen in Fig. [5J A more elabo 



rated version of this figure can be found in lDecin et al 
( 2000h . for G-type stars and cooler, differences in stel- 



lar parameters within one spectral sub-class (e.g. from 
K1III to K2III) may cause differences in the (broad) 
NIR molecular absorption bands in the order of 5 - 8 % 
(at a resolution of 100), an uncertainty that should be 
directly translated to the uncertainties attributed to 
the template spectra. For hotter stars, the deviations 
are smaller and arc mainly seen in the strength of the 
hydrogen lines. The continuum flux may differ by up 
to -4%. 

Recently, lEngelke et al. found, for a subset of 

33 standards with spectral types between G2 and M5, 
available in the CWW network, that the derived spec- 
tra (by Engelke and co-workers) display systematic dif- 
ferences with the composites/templates in the CWW 
network, in that Engelke standard spectra have 4-7 % 
lower flux in the 1-4 spectral region. The spectra of 
Engelke et all (|2006h remove the A-K star calibration 



bias recently noted in the calibration of Spitzcr/IRAC, 
which was based o n templates in the CWW network 
( Reach et al1l2005l) . 
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Table 1. Summary of uncertainties attributed to the standard star spectra in the CWW network in the 2-200 /im 
range. 



Description 



Uncertainty Spectral Type Wavelength Region 



propagation of systematic errors 
e.g. due to use of Vega as prime calibrator 
— > absolute flux calibration of Sirius w.r.t. Vega 

use of Engelke function (Engelkc 1992) 

creation of 'template' spectrum 
— > molecular features 
— > continuum 

presence of chromosphere/ionised wind 
presence of circumstellar dust 
systematic difference with new Engelke 
standards (Engelke et al. 2006 ) 



1-3% 


A-M 


2-200 urn 


3% 


A-dwarfs 


22-35 /im 


8% 


K-giants 


22-35 /im 


up to 8 % 


G-K 


around 2.3, 4.0, 4.2, 8/im 


up to 4 % 


A-M 


2-200 /im 


£10% 


G-M 


A > 100 /im 


£10% 


A-M 


A > 2 |Um 


4-7% 


G2-M5 


1-4 /im 



The net result is that the CWW network has pro- 
vided the astronomical community with a consistent net- 
work containing some 600 stars. Calibration experiments 
for the SPIRIT III infrared instrument onboard the MSX 
have demonstrated that in the 6 (broad) spectral bands 
(with wavelength between 4.29 and 21.34 /im) the absolute 
flux calibra tion for 8 compo sites in the CWW network 
is < 2.3% (jPrice et al.l 120041 1. However, previous discus- 
sion also shows that for some spectral regions of the IR 
SED, the uncertainties attributed to the template spectra 
or extrapolated continuum spectra in the CWW network 
should be larger than now quoted (for an overview, we re- 
fer to Table Q]) . A typical example of the attributed total 
uncertainty for two template spectra is given in Fig. [5l In 
the spectral regions where molecules absorb photons sub- 
stantially, at least an 8 % uncertainty (see above) should 
be attributed to the tcmplatcd spectra. Another exam- 
ple is the (unknown) contribution of a (possible) chro- 
mosphere, ionised wind or dust-driven wind, which may 
yield un certainties in the order o f 10% at 100 /zm (see 
Sect. 14.31 van der Bliek et al. 19961 ). much higher than the 
5.67% uncert ainty attribut e d to th e extrapolated compos- 
ite spectra bv lCohen et al. ( 1996b ). 

In order to remove the 'molecular feature' bias that can 
enter into the template construction, the computation of 
ab-initio theoretical model atmosphere spectra may offer 
a solution (see next section). Moreover, theoretical atmo- 
sphere spectra can be computed at any spectral resolution 
(comparable to that of the instrument) and in that way 
compete with the template spectra, which have quite low 
spectral resolution. 

3.4. Model atmosphere spectra 

Different groups have given a lot of effort into develop- 
ing trustworthy hydrostatic 1- dimensional (1-D) model 
atmosphere codes, the most famous ones in the field 
of cool stellar atmospheres (T e g ^ 10 000 K) being 
the ATLAS atmosphere code (jKurucz ll970L Il993l Il996l 
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Figure 5. Upper panels: Template spectra of £ Dra 
(K2 III) and a Tuc (K3 III). The composite spectrum of 
a Boo (K2 IIIp) has been used to construct the template 
of £ Dra (K2 III), and for a Tuc (K3 III) the composite of 
a Hya (K2 II-III) is used. Lower panels: the total uncer- 
tainty attributed to the templated spectra bv lCohen et al 
(1999). Vlain absorption features are the CO fundamen- 



tal around 4.5 /im, CO first overtone around 2.3 /im, SiO 
fundamental around 8 /im and SiO first overtone around 
4.2 /im. 



Castelli fc Kurucd 2004 ), the MA RCS code (originally pre- 



sented in lGustafsson et al.lll975l. with t he current grid bo 



ing d escribed inlGustafsson et al.ll2003 ). and the phoenix 



code ( Hauschildt et alj 1999al lb|) Also in the case of late- 
type (variable) asymptotic giant branch (AGB) stars, the 
1-D hydrostatic models are still useful in providing us 
with physical information, due to the great detail now 
included in their treatment of opacities. Hydrodynamic 
models in 1-D (e.g iBessell et al.|[l996l IWinters et~al1l2000l: 
H5fner et alj 120031 ) include detailed physics, like time- 
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dependent dust formation and the coupling between gas, 
dust, and radiation. Few dynamical models are now able 
to simultaneously solve the equations of hydrodynam- 
cics and frequency-dependent radiative transfer leading to 
cons istent dynamical de nsity-temperature structures (see, 
e.g- lHofner et alJl2003h . Us ins these kinds of models, one 
can already quantitatively reproduce the variation in line 
profiles due to the influence of gas velocities (see, e.g., 
Nowotnv et al. I l2005h . Recently, 3-D hydrodynamical sim- 
ulations of stellar surface convection have become feasible 
thanks to advances in computer technol ogy and efficient 
numerical algor ithms ((Collet et al.ll2006l ). The 3-D mod- 
els can shed light on the coupling between convection and 
pulsation, as well as be employed in element abundance 
analysis. These 3-D model computations are, however, still 
very computer-time intensive, and it is not known how 
well they represent the surface layers of the stars, so spec- 
trum synthesis of a whole sample of stellar spectra is still 
beyond the scope. 

For the scope of this project of developing accurate 
SEDs for IR stellar standards, we therefore rely on a 1-D 
hydrostatic model atmosphere code. We here opt to use 
the MARCS code (see next section). 



and phoenix model fluxes. Moreove r ^ the MARCS mod- 
els are based on the lAsplund et al.l ((200J) solar abun- 
dances, while the ph oenix models follow the results of 
Grevesse et al.l (|l996T ). Note also that the marcs model 
uses a more complete OH and SiO line list at A > 10 fim 
(see also Fig. [7]). 
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4. Status and accuracy of marcs model 
atmosphere spectra in the IR 

The MARCS model atmosphere and spectrum synthe sis 
code was developed in Uppsala ( Gustafsson et al. 2003 [H 
This local thermodynamic equilibrium (LTE) model at- 
mosphere code is built on the assumptions of spherical 
or plane-parallel stratification in homogeneous stationary 
layers and hydrostatic equilibrium. Energy conservation 
is required for radiative and convective flux, where the 
energy transport due to convection was treated through a 
local mixing- length theory. For a discussion of the method 
for solving the radiative transfer equations in the at- 
mospheric mode ls a nd spectrum synt hesis, we refer to 
Nordlundl <|l984h and lPlez et all ((l992h . 

In the framework of the release of a grid of MARCS 
models, colleagues from the University of Uppsala are 
comparing the model structures of the MARCS, ATLAS, 
and phoenix codes. An example typical of a cool K gi- 
ant is shown in Fig. [6l The plot is made in the \ 2 F y - 
space to clearly distinguish the differences. The contin- 
uum level between the marcs and phoenix models do 
agree within 0.5%. A larger difference of ^5% is seen 
between the MARCS and ATLAS models for the contin- 
uum at A > 10 /im. The difference in the molecular fea- 
tures is 554% between MARCS and phoenix and <8% 
between MARCS and ATLAS at a resolution of 100. The 
reason for the difference is the use of a plane-parallel 
geometry in the case of the ATLAS models, whereas the 
phoenix and MARCS models use a spherically symmetric 
geometry. Slightly different surface temperature and geo- 
metrical extension explain the difference between MARCS 



Figure 6. Upper panel: Comparison between the ATLAS, 
phoenix, and MARCS model spectra. The PHOENIX model 
is shifted downwards with a factor 10, and the MARCS 
model with a factor 100. Middle panel: Comparison be- 
tween the ATLAS, phoenix and MARCS model spectra at 
a resolution of 100 for A < 10 /im. At longer wavelengths, 
only the 8 wavelength points as given in the ATLAS models 
(at 20, 40, 60, 80, 100, 120, 140, and 160 fim) are retained. 
Difference between phoenix and MARCS model fluxes is 
barely visible. Bottom panel: Ratio between MARCS and 
ATLAS (black) and phoenix (gray) model fluxes. 



For spectrum synthesis, data on the absorption by 
atomic and molecular species are collected from dif- 
ferent databases. For a discussion of various available 
atomi c and molecular line lists, we refer to 
( 2000l ). For the purpose of this work, the 



Decin 



follow- 
ing infrared spectrosc opic line lists were us e d: CO 
line list computed by Goorvitch fc Chackerianl ( 1994T ) . 
SiO by lLanghoff fc Bauschlicherl (Il993h. CN bv Plez 
(priv. comm.), OH by Goldman et al. ( 19981 ). H 2 by 
Partridge fc Schwenkef ((l997|), NO, HF, NH, HC1, and 
C H by Sau v al (p riv. comm.), and atomic line lists 



1 http://marcs.astro.uu.se 



by Ivan Hooj (|l998h . by iHirata fc Horaguchil (Il995l). by 
Sauval (priv. comm.), an d of VALD (IPiskunov et al.lll995 : 
iRvabchikova et al 1997 : Kupka et alj 19991 ). Using the 
model photosphere as input, synthetic spectra are calcu- 
lated for a typical resolution of A/AA ~ 300 000, even 
though the final instrumental resolution is often lower. 
With a typical microturbulcnce of 2 km s -1 , this means we 
are certain to sample all lines in the atomic and molecu- 
lar database. We note that, for the purpose of calculating 
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theoretical spectra in the mid to far-IR, some line lists 
are still far from complete or accurate (see Sect. I4.5|) . An 
typical example of a K2 giant is shown in Fig. [71 At a 
medium resolution of A/AA = 1500, the depression of flux 
due to the line veiling is ~3 % at wavelengths longer than 
~30 fim. 

For the purpose of using model atmosphere spectra to 
represent the SED of standard stars, an appropriate anal- 
ysis of different sources of uncertainties contributing to 
total uncertainty in the spectrum predictions is in place. 
The following sections discuss the effects of sources of er- 
ror as a) the input stellar parameters, b) uncertainties in 
the model assumptions, c) the possible presence of a chro- 
mosphere, ionised wind or circumstcllar dust shell, d) the 
adopted continuum opacity, and e) the used line lists. 



4.1. Dependency on stellar parameters 



Model atmospheres arc defined by the fundamental pa- 
rameters effective temperature, gravity, and metallicity 
(and stellar mass or radiu s in case of a sphe rical geom- 
etry). As demonstrated in iDecin et al.1 ( 200fj| ). the influ- 
ence of the stellar mass on the synthetic spectrum is small. 
Varying the other fundamental parameters shows the ef- 
fect of errors in the determination of the fundamental pa- 
rameters on the synthetic flux distribution. Although this 
effect is dependent on the full set of fundamental stellar 
parameters, we may summarise that, for A > 50 fim (1.), 
varying the effective temperature by ~200K for stars with 
spectral type between G and K, roughly corresponds to a 
change of ^4% in the continuum flux in the IR, (2.) a 
change in the logarithm of the surface gravity of 0.20 dex 
introduces uncertainties in the continuous flux distribu- 
tion of about 0.5 %, and (3.) an uncertainty in the metal- 
licity of about 0.20 dex corresponds to an uncertainty in 
the IR continuum flux of ~0.1 %. The uncertainty in the 
near-IR and on molecular absorption features is, however, 
much greater, as can easily be seen in Fig. [8j arising to 
7% (at a resolution of 100). 



The same order of magnitude on the uncertainties also 
applies to the continuum predictions of the hotter A-type 
stars for the same relative changes in stellar parameters. 
On one hand, one may argue that the absence of broad 
molecular absorption bands makes these hotter dwarfs 
more eligible as standards. On the other hand, exactly 
the presence of molecular features in the cooler giants can 
be used as a strong diagnostic tool to estimate the stellar 
parameters, and hence to provide a fiducial prediction of 
the whole IR spectrum. Especially in spectral regions with 
strong hydrogen lines, for which the computation of the 
self-broadening rem ains problematic ( Barklem et al.ll2000t 
Decin et al.ll2003b ). the use of cooler giant as stellar stan- 
dards will improve the accuracy of the spectrophotomctric 
calibration. 




[Fe/H]=-0.50 / [Fe/H] = -0.30 J 



10 100 
wavelength [jum] 

Figure 8. Ratio between synthetic spectra with differing 
stellar parameters at a resolution of A/AA = 100. The 
standard model has the same stellar parameters as in 
Fig. UJ Upper panel: Ratio between a synthetic spectrum 
with T efI = 4320K and with T cff = 4520K. Middle panel: 
Ratio between log g= 1.50 and log g= 1.70. Bottom panel: 
Ratio between [Fe/H] = -0.50 and [Fe/H] = -0.30. 



4.2. The temperature distribution T(t) 

Model atmosphere calculations are based on a number of 
assumptions, one of them being radiative and convective 
flux conservation also in the outermost layers of the pho- 
tosphere. This simplification may be the reason that the, 
otherwise almost perfect, match between the FTS-Kitt 
Peak Spectrum of a Boo and theoretical predi ctions based 



on MA RCS atmospheres (see Figs. 16-17 in IDecin et al 



2003bl ) shows small deviations at the 1-2 % level at a res- 
olution of A/AA ~ 60 000 in the low-excitation CO and 
OH lines. 

Not only for K giants (as a Boo), but also for the 
most studied star, the Sun, there is still quite some de- 
bate concerning the assumed mean thermal profile in the 
outer layers of the photosphere. Both for G and K-type 
stars, indications are found from, e.g., the Ca II H and K 
lines that the temperature structure has a minimum before 
segueing into the mec hanically heated chromosphere (e.g. 
Avres fc Linskvlll975l ) . Controversially, the analysis of CO 
lines indicate s a cooler brightness temperature at the same 
altitude fe.g.lAvres fc Testermanlll981uWiedemann et al 



1994 lAvres et all 12006}) . 



As done by Ivan der Bliek et al 



(Il996h and 

Vanhollebekej (|2003l ). one can simulate uncertainties 
in the temperature structure T(t) to mimic different 
effects: a) flatten or steepen T(r) for a Rosseland optical 
depth r ross < 0.01 to study the effect of less, respec- 
tively more, line blanketing, b) steepen T(r) around 
r ross = 0.1 to study the effect of convective overshoot in 
the outer layers of the model photosphere, or c) steepen 
T(t) for t ioss > 1 to study the effect of convection in 
the deeper layers of the atmosphere. From Table 3 in 
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Figure 7. Contribution of CO (blue), 
SiO (orange), OH (green), CN (pink), 
H2O (red), atoms (purple) to the 
2 — 200 /urn synthetic spectrum cal- 
culated from the atmosphere model 
with stellar parameters T c ff = 4320 K, 
log g [cm/s 2 ] = 1.50, mass M = 
1.1 M Q , [Fc/H] = -0.50, microturbu- 
lence & = 2 km s" 1 , e(C) = 7.96, 
e(N) = 7.61, e(0)=8.68, e(Mg) = 7.33, 
e(Si) = 7.20, and 12 C/ 13 C = 7. The first 
six panels display the contribution at a 
resolution AA = 0.5 A in color and at 
a medium resolution A/ A A = 1500 in 
black. The total contribution to the full 
spectrum is shown in the seventh panel 
at the high resolution of AA = 0.5 A, 
and in the bottom panel at the medium 
resolution of A/AA = 1500. Note the 
lack of atomic absorption features for 
wavelengths longer than ~50 /im. 



van der Bliek et al. (1996) and the results of lVanholleb ckc 
( 20031 ). we can conclude that this type of error in the 
temperature distribution gives rise to uncertainties in the 
predicted continuum flux at 100 /im lower than 3.5%. 
This uncertainty can, however, be strongly reduced in 
case high-resolution spectra arc available for atomic or 



molecular lines with different strengths and excitation 
energies. From a proper analysis of the line shapes, one 
can pin down T(r) to ~50K in the outer layers, reducing 
the unc ertainties in IR continuu m flux predictions to 
-1-2% (Ivan der Bliek et al.|[l996T). 
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4.3. The presence of a chromosphere, ionised wind, or 
circumstellar dust shell 

One of the largest uncertainties in model spectrum 
predictions is the (unknown) presence of a chromo- 
sphere, an ionised wind, or a circumstellar dust en- 
velope. These structures may easily yield a flux ex- 
cess in t he order of 10% a t mid to far-infrared wave- 
lengths (Ivan der Bliek et al. 



, 19961) . Recent studies by 
Dehaes et al.l ( 20071) show that all of the six studied cool 



standard stars, used for the calibration of the Short- 
Wavelength Spectrometer (SWS, 2.38 - 45.2 p,m) onboard 
ISO, with spectral types between K2 III and MO III show 
a flux excess due to an ionized wind. Only for two of them, 
the flux excess starts at wavelengths shorter than 1 mm, 
but still longer than ~200 fxm. In the case of the Spitzer- 
IRS (5.3 - 18 /im) calibration, the IRS team had to reject a 
considerable fraction of candidate standards as calibrators 
due to IR excess emission (Sloan 2005 , priv. comm.). 

4.4. The continuous opacity source: 

The continuous opacity of solar-type and cooler stars is 
dominated by free-free absorption of the negative hydro- 
gen ion. The MARCS model atmosphere calculations make 
use of the free-free a b sorpt ion coefficients as calculated 
by iBell fe Berrington (1987). A quantitative asses sment 
of the reliability of absorption coefficients by Ijohn 
(Il994l) shows that th e abso rption coefficients tabulated 
by IBell fc Berrington! (jl987l ) are accurate to about 1 % 
for wavelengths greater than 0.5 /xm over the tempera- 
ture r ange between 1400 and 10 080K. IStillev fc Callaway 
( 1970f ) neglected the adiabatic exchange potential term 
yielding absorption coefficients that may be off by ~4.5% 
at temperatures lower than ~2500K. Using these coeffi- 
cients, one may introduce errors in the order of 4 % at 
wavelengths < 10 /im, diminishing to less than 0.3% for 
longer wavelengths (see Fig. The prominent feature, 
seen around 8 /im in Fig. [5J arises from the response of 
SiO on the slightly different temperature structure. 

4.5. Line lists 

Atomic and molecular data bases (including very differ- 
ent numbers of molecular lines and based on computations 
spanning a wide range in quality) are used in the construc- 
tion of model atmospheres. Knowing the limitations on the 
accuracy and completeness of the used line list is of key im- 
portance when calculating theoretical spectra. The MARCS 
code, with the input data used, is particularly tuned for 
computations of late-type stars where molecular opaci- 
ties play an important role. As said in the i ntrodu c tion o f 
Sect. |H we here rely on the study done by iDecin ( 2000h . 
In this study, it has been shown that for molecules such 
as CO, SiO, and OH, different molecular line lists reach 
a high level of accuracy and completeness. In the case of 
CN, the dissocia tion energy is, however, still a matter of 
debate (see, e.g.. lLambertlll994l ). Most problematic is the 
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Figure 9. Ratio between the synthetic spectrum 
calculated with the H g absorption coefficients of 



Bell fc Berrington] ( 1987 ) and the spectrum based on th e 



H ff absorption coefficients of Stilley fc Callaway ( 197Clh 



The same stellar parameters as in Fig. [7] are used, data 
are reduced to a resolution of A/AA = 100. 

situation for the H?0 line lists. The most often used theo- 



retical H2O data bases are the SCAN list (lJ0rgensen et al. 
200lh , the list of IPartridge fc Schwenkel .19971). and the 



very recent list of iBarber et al.l (|2006) . These line lists 
are likely to be satisfactory for opacity calculations and 
probably for comparing observed and synthetic spectra at 
low and medium resolution, but at higher resolution line 
positions based on la boratory measurements (e.g., from 
Polvanskv et al. 19971 ) should be used if possible. 

Focusing on the wavelength range between 60 and 
210 ^tm (covered by the ESA PACS instrument onboard 
Hcrschel, launch foreseen in 2008), one should realise that 
the atomic VALE0 database presently only contains 1380 
lines with A < 122 /im, the NIST d atabase^ o n ly 90 atomic 
lines, while the atomic line list of van Hoo3 ( 19981 ) tabu- 
lates 13 527 lines, of which the oscillator strengths are only 
known for 1660. 

MARCS mode l atmo spher es were used by, among o th- 



Decin et al.l (|2000l) and iRvde fc Eriksson! (|2002[) to 



ers, 

predict IR spectra of cool stars. In their comparison 
between the high-resolution (A/AA ~ 60 000) Fourier 
Transform Spectrometer (FTS ) spectrum o f a Bo o and 
theoretical MARCS predictions. IDecin et al.l ( 2003bJ ) show 
that, in the case of this well-known giant, the model pre- 
dictions in the 0.9 - 5.3 /im range only differ from the 
observational data by 1-2 %! In general, IR model spec- 
tra predictions in the wavelength range < 25 /im for A-G 
dwarfs and for non-pulsating K-M0 giants may be as ac- 
curate as 3 % at medium resolution (A/AA ~ 1500) and as 
~5 % at high resolution (A/AA ~ 1 00 000) , as judged fr om 
the available observational data ( Decin et al.l l2003bl lal). 



2 http:/ /ams. astro. univie.ac.at/vald/ 

3 http:/ /physics. nist.gov/cgi-bin/AtData/main_asd 
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An exception may be the hydrogen-line predictions due 
to the problematic c omputation of the self-broadening 
( Barklem et al. l2000h . At longer wavelengths, the accu- 
racy and resolution of today's modern instruments remain 
too poor to constrain the model atmosphere spectra at a 
few percent level. Luckily, the depression of flux due to the 
line veiling in the FIR is estimated to be ^ 3 % at a reso- 
lution of 1500 (sec Fig. [7]), rendering the representation of 
the SED of stellar standards by theoretical model atmo- 
sphere spectra still useful. In order to avoid the propaga- 
tion of inaccurate predictions of molecular/ atomic lines to 
the RSRF determination, different spectral classes should 
be used. At temperatures lower than 3500 K, H2O becomes 
a dominant opacity source, excluding giant stars with a 
cooler spectral type than MO III as standard calibrators. 

4.6. Summary 



(GBPPfl as obtained bv lHammerslev et alJ (| 19981 ) and 
Hammerslev &; Jourdain de Muizorl "( 2003 ). 



For the MARCS atmosphere s pectra, stellar parameters 



as determined by iDecin et al. ( 2003a ) were used as in- 
put parameters. Main parameters for the absolute flux, 
T e fj, and 6d, are listed in Table [3] For the computation 
of the flux values of the Engelkc function for both com- 
posites, T ff and 6d, as given in the headers by CWW, are 
used. Note that for (3 And, the re-scaled (by CWW) angu- 
lar diameter of 13.71 mas is used: based on t he InfraRed 
Flux Method (IRFM). lBTackwell et all (|l99lh determined 
a lower angular diameter of 13.219 mas corresponding to 
the stellar temperature of 3839 K. As noted in Sect. 13.21 
the Engelkc function needs a higher angular diameter (or 
T ff) for low-gravity, low-temperature giants in order to 
attain the correct flux level. For both templates, same val- 
ues as for the MARCS spectrum are used. The input for 
the blackbody calculation is the same as for the Engelkc 
function. 



In last few sub-sections many uncertainties on the 
calculations of model atmospheres are summarised (see 
Table [2] for an overview), and the total uncertainty is not 
(as often done) estimated as the root-sum-square of all 
above-mentioned uncertainties since they are often mutu- 
ally dependent. The total uncertainty on the stellar refer- 
ence SEDs depends on the spectral type of the target, 
the wavelength range under study and the instrumen- 
tal resolution. In case the standard star has been prop- 
erly studied using high- resolution optical and/or near-IR 
data and the presence of an IR excess can be excluded, 
one can constrain T e ff and the temperature distribution 
to within ~50K, yielding near-IR (continuum + line) 
medium-resolution predictions better than 1-2 %, mid-IR 
predictions better than ~3 % and far-IR predictions bet- 
ter than ~5 % for stars with earlier spectral type than M0. 
Each of the discussed spectral classes (early A dwarfs, so- 
lar analogs, and G9-M0 giants) has its own drawbacks. 
Consequently, with the aim of constraining the RSRF at 
a few percent level, the selected stellar standard candles 
should cover a wide diversity in spectral types to avoid 
biases to spectral features typical of different spectral 
classes. 

5. Comparison between different IR stellar 
reference SEDs 

It is instructive to compare the different IR stellar ref- 
erence SEDs presented in this paper. We therefore have 
chosen to compare the proper blackbody, Engelkc func- 
tion, and model atmosphere spectrum with the spectra 
of two composites (as explained in Sect. 13. 3[) (a Boo 
(K2 IIIp) and /3 And (M0 III)) and two templates (£ Dra 
(K2 III) and a Tuc (K3 III)) from the sample of CWW 
(see Fig. [10]). Figure [10] also plots the IRAS LRS and 
PSC data and different photometric data points collected 
from the ISO Ground Based Preparatory Programme 



Table 3. Stellar parameters used to for the calculation 
of (1.) the MARCS model atmosphere spectra (columns 
2-5) and (2.) the Engelke function and blackbody (BB) 
(columns 6-7) in Fig. [TO] Effective temperature, T e ff, is 
given in Kelvin, gravity, g, in cm/s 2 , and angular diameter, 
6d, in mas. 



target 


T c ff 


MARCS 
logg [Fe/H] 


e d 


Engelke and BB 
T c ft 9d 


£ Dra 


4440 


2.40 


0.10 


3.09 


4440 


3.09 


a Boo 


4320 


1.50 


-0.50 


20.66 


4362 


21.12 


a Tuc 


4300 


1.35 


0.00 


5.90 


4300 


5.90 


(5 And 


3880 


0.95 


0.00 


13.48 


3839 


13.71 



Concentrating on the two composites, a Boo and (3 
And, of which a large part of the CWW spectrum con- 
sists of observational data, it is immediately clear that the 
Engelke function does a much better job than the black- 
body function of representing the observed continuum 
stellar SED. The simple blackbody function can not even 
be used to estimate broad-band features in the RSRF. 
Between ~23 and ~ 35 /im, the CWW spectrum of both 
composites is given by their Engelke function. As discussed 
in Figs.Q]and[2] the shape of the continuum flux differs be- 
tween a proper MARCS model spectrum and the flux values 
given by the Engelkc function, the reason being the neglect 
of the influence of the gravity and sphericity effects, which 
are important for giant stellar atmospheres. Between 20 
and 200 /im, the maximum difference between the Engelke 
function used by CWW and the MARCS model atmosphere 
for a Boo is near 60 /im, where it attains ~2.5%, rising 

4 http:/ /www.iso.vilspa.esa.es/users/expLlib/ISO/wwwcal/ 
isoprep /gbpp /photom / 
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Table 2. Summary of uncertainties attributed to the theoretical atmosphere spectra in the 2 200 /mi 
range at a resolution A/AA ~ 100. 



Description 


Uncertainty 


Spectral Type Wavelength Region 


• dependency on stellar parameters 








— > molecular features 


up to 8 % 


G-K 


around 2.3, 4.0, 4.2, 8 /mi 


— > continuum 


up to 4 % 


A-M 


2-200 /tm 


• uncertainties on T(r) 








continuum flux (without high-resolution data constraints) 


<>3.5% 


A-M 


2-200 /im 


continuum flux (with high-resolution data constraints) 


1-2% 


A-M 


2-200 /im 


• presence of chromosphere/ionized wind 


£10% 


G-M 


A > 100 /im 


presence of circumstellar dust 


£10% 


A-M 


A > 2 /mi 


• continuous opacity by 


1% 


A-M 


2-200 /im 


• line lists 


<,3% 


A0-M0 


2-200 /im 


OVERALL budget: for approved standards 
with high-resolution data constraints 


1-2% 
-3% 
-5% 


A0-M0 
A0-M0 
A0-M0 


near-IR 
mid-IR 
far-IR 



to ~7.5 % when the same T e g and 9d are used as for the 
MARCS model spectrum (Fig. [TTj) . 

The large absolute difference between the IRAS LRS 
spectrum of a Boo and both the MARCS model atmo- 
sphere spectrum and CW W composite (~ 35 %, sec also 
Van Malderen et al.ll2004l ) remains unexplained. The LRS 
raw data are extracted from the Groningen IRAS database 
and calibrated with the lrscal routine in the GIPSY 
package. Although the IRAS LRS data were originally 
not meant to be absolutcd ca librated, absolute cal ibra- 
tion factors were det ermined bv lVolk fc Cohe n (1989!) and 



Cohen et al 



Cohen et al 




The same procedure was used by 
who however only needed a factor 
of 0.95 to splice the LRS data to the CGS3 data between 
7.5 and 13 /im. 

Inspecting the molecular absorption features in both 
templates, £ Dra and a Tuc, a clear difference is visible be- 
tween the template spectrum (red) an d the MARCS mode l 
atmospheres (black). As discussed in Decin et al] ( 2000h 
and in Sect 14.11 differences in abundance pattern (here 
mainly the C and O abundance), temperature, gravity, 
and metallicity result in substantial differences in strength 
in molecular absorption. In the case of £ Dra, the compos- 
ite of a Boo was used to construct the template, the com- 
posite of a Hya (K2 II-III) was used for a Tuc. Inspecting 
the literature study of a Boo an d £ Dra presented in 
Appendix D in lDecin et al.l ( 2003af ) clearly shows that in- 
dependent studies incorporating both objects in general 
yield an effective temperature higher by ~120K, a loga- 
rithm of the gravity higher by ~ldex, and a metallicity 
[Fe/H] higher by ~0.50dex for £ Dra with respect to a 
Boo. Figure [12] compares the MARCS flux predictions of 
a Boo and £ Dra. Especially in the regions of molecular 
absorption, the (relative) difference between both spectra 
may be 10 % or even larger, explaining why the difference 
between the CWW spectrum and the MARCS model spec- 
trum is larger than the uncertainties quoted by CWW. 




MARCS model Alpha Boo 
Engelke function (T eff =4362. S d =21.12 mas) 
Engelke fjnc;ion (T eff =4320, fl d =20.72 mas) 




wavelength [yum] 

Figure 11. Upper panel: Comparison between the MARCS 
model spectrum of a Boo (A/AA = 100) and predictions 
using the Engelke function with parameters T e ff = 4362 K 
and 9d = 21. 12 mas (dashed line), and T e ff = 4320K and 
9d = 20.72 mas (dotted line). Bottom panel: Ratio between 
MARCS model predictions and both Engelke functions. 



Consequently, the use of the composite spectrum of the 
metal-deficient K2-gaint a Boo template to repre- 
sent a whole class of K2 giants in the CWW network will 
introduce additional uncertainties, which may propagate 
through the RSRF determination. 

For bot h composites , the fa r-IR extrapolations con- 
structed by Cohen et al. ( 1996bl ) differ by a few percent 
from the MARCS atmospheric predictions, since the MARCS 
models are computed in spherical geometry, while the 
CWW predictions are interpolations in a grid of plane- 
parallel models. For the cooler, more-extended giant cali- 
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Synthetic spectrum 

Cohen SED 

Cohen SED +/- a 

Cohen photometry 
IRAS LRS 



IRAS PSC (colour corr.) +/- a 



Engelke Function 
Blackbody 



Figure 10. Comparison between the 
various reference SEDs discussed in 
this paper in the 2 to 200 fj,m wave- 
length range, a Boo (K2 IIIp) and (5 
And (MO III) represent two compos- 
ites from the CWW network, and £ 
Dra (K2 III) and a Tuc (K3 III) rep- 
resent two templates. MARCS model at- 
mosphere spectra are displayed at a res- 
olution A/ A A = 100. IRAS LRS and 
PSC data and several photometric data 
points collected from the ISO GBPP 
are also shown. 



brators, sphericity effects may yield an infrared flux excess 
in the order of ~5 % (see Fig. [2|. 



6. Theoretical atmosphere spectra for stellar 
candles and how to obtain them 

The four reference SEDs constructed from the MARCS at- 
mosphere code presented in Fig. [10] are available via the 
electronic version of the article. Other reference SEDs 
constructed in the framework of the calibration of e.g. 
ISO, Spitzer, MIRI etc. are available upon request. If one 
needs extra reference SEDs for additional stellar calibra- 
tion sources, one may always contact the authors, who 
are willing to provide you with the appropriate theoreti- 
cal model atmosphere spectrum in the wavelength range 
and at the resolution requested. 



7. Conclusions 

In this paper, various sets of reference SEDs used for deter- 
mining the spectrophotomctric calibration of IR spectrom- 
eters (onboard satellites) have been discussed. Our main 
emphasis was on the stellar reference SEDs, with special 
focus on model atmosphere spectra. It is shown that the 
predicted medium-resolution IR model atmosphere spec- 
tra are accurate within 1 -2 % in the near-IR, ^3 % in the 
mid-IR, and ^5 % in the far-IR for stars with spectral 
types earlier than M0. From the four types of stellar ref- 
erence SEDs discussed in this study (blackbody, Engelke 
function, templates in the CWW network, and theoretical 
atmosphere spectra), it is believed that theoretical atmo- 
sphere spectra make up the best representations nowa- 
days for the stellar SED, especially in case one wants to 
calibrate instruments with a spectral resolution ^ 500. 
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A [/im] 

Figure 12. Upper panel: Comparison between the MARCS 
model spectrum of a Boo and £ Dra at a resolution of 
1500. Both spectra have a fictitious angular diameter of 
1 mas. Bottom panel: Same as for the upper panel, but 
now at a resolution of 100. In regions of strong molecular 
absorption, the maximum difference is 15% for A/AA = 
1500, and 8% for A/AA = 100. 



We note, however, that at A > 20 /mi, stellar calibrators 
might be too faint to be used in the spectrophotometric 
calibration pedigree. 

Since the ultimate goal of the calibration system is 
often to be capable of deriving spectral flux values that 
are trustworthy to 3 % or better on both absolute and 
relative scales, one should aim at building a highly accu- 
rate system of stellar reference SEDs. A good sky cover- 
age by the calibrators is an important ingredient in terms 
of a time-efficient determination of the (spectrophotomet- 
ric) calibration. In that sense, the IR network constructed 
by CWW is a good starting point. However, the CWW 
network has its limitations in terms of accuracy at rep- 
resenting the (molecular) spectral features and the SED 
at wavelengths longer than ~20 ^m. For a list of can- 
didate calibrators, one therefore should put effort into 
obtaining ancillary observational data to both (1) con- 
strain the reliability of the candidates as standard stel- 
lar sources and (2) estimate the stellar parameters to 
compute a set of highly reliable theoretical model at- 
mosphere SEDs. This set should compromise standard 
stellar candles with different spectral types ranging be- 
tween AO V and M0 III. Experiences with the Short- 
Wavelength Spectrometer (SWS) onboard ISO and the 
InfraRcd Spectrometer (IRS) onboard Spitzer have shown 
that this set should consist of some ten to fifteen stellar 
calibrators. In the framework of the calibration plans for 
PACS onboard the ESA-Herschel satellite and MIRI on- 
board the NASA-JWST satellite, a set compromising the- 
oretical spectra of ~15 stellar calibrators is in preparation 
(jDecin et al.ll2OO70 . The input stellar parameters for each 
set of model spectra have to be estimated from a proper 



analysis of high-resolution optical or near-IR data. A sys- 
tematic cross-calibration with planets and asteroids is the 
last step in the spectrophotometric calibration pedigree. 
Only then can a highly accurate calibration system be de- 
veloped. 
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